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SECTION  I 


INTRODUCTION 

The  concept  of  eddy  viscosity  has  been  applied  to  predict  turbulent 

fl  2  31 

boundary  layers  in  hypersonic  flow  regions.  *  '  Hopkins  et  al.  found 

(41 

that  the  multilayer  eddy  viscosity  model  by  Cebeci'  '  underpredicted  the  skin 
friction  coefficient  by  around  10f  for  hypersonic  flow  conditions.  Extension 
of  the  essentially  incompressible  eddy  viscosity  model  to  compressible  flow  is 
based  on  the  fact  that  the  turbulence  producing  mechanism  remains  similar  in 
spite  of  the  presence  of  temperature  or  density  fluctuations.  The  prediction 
of  compressible  turbulent  boundary  layers  by  eddy  viscosity  models  indeed 
yields  essential  agreement  with  experimental  data  up  to  a  Mach  number  of 
five. However,  in  the  hypersonic  flow  region  Bushnell  and  Beckwith 
revealed  a  dependence  of  the  Reynolds  stress  on  the  density  fluctuations.  In 
fact,  Bushnell ^  introduced  a  mixing  length  model  for  the  density  fluctuations 
in  his  formulation  of  the  apparent  viscosity  coefficient.  A  few  specific 
comparisons  were  made  to  illustrate  the  effect  of  the  fluctuating  property  on 
the  mean  flow. 

Recent  experimental  data  recorded^® at  high  Mach  numbers  reveals 
significant  pressure  fluctuations  in  the  turbulent  boundary  layer.  A  common 
characteristic  of  the  measured  pressure  data^’^  indicates  a  distinct  nonral 
pressure  gradient  within  the  turbulent  boundary  layer.  At  the  present  time, 
no  analytical  method  has  been  developed  to  include  this  phenomena. 

Numerical  solutions  of  the  compressible  turbulent  boundary  layer  with 
heat  transfer  are  commonly  obtained  by  assuming  a  constant  turbulent  Frandtl 
number.  The  Prandtl  number  is  generally  assigned  to  be  near  unity. 

Recently,  extensive  experiments  have  been  devoted  to  studying 
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the  Prandtl  number  distribution  across  the  turbulent  boundary  layer.  The  work 

19) 

on  the  turbulent  Prandtl  number  by  Meier  and  Rotta'  '  revealed  that  the 
turbulent  transport  of  energy  decreases  more  rapidly  toward  the  wall  than  the 
momentum  transport.  Therefore,  the  turbulent  Prandtl  number  exceeds  the  value 
of  unity  near  the  wall  but  decreases  slightly  within  the  sublayer  region. 

Their  findings  fall  within  the  uncertainty  which  envelopes  the  turbulent 
Prandtl  number  distribution  by  Simpson  et  al.^^  Their  results  seem  to 
indicate  that  the  heat  transfer  and  temperature  profile  of  turbulent  flows 
cannot  be  ascertained  by  an  oversimplified  constant  turbulent  Prandtl  number. 
Again,  no  systematic  evaluation  of  the  effect  of  variable  Prandtl  number  on 
the  heat  transfer  and  temperature  distributions  in  a  boundary  layer  has  been 
performed. 

The  present  analysis  Intends  to  Investigate  the  hypersonic  turbulent 
boundary  layer  with  heat  transfer.  The  explicit  dependence  of  the  density 
fluctuation  Is  Introduced  Into  the  eddy  viscosity  model  in  a  similar  manner 
to  that  of  Bushnell  et  al.,  except  that  the  density  fluctuation  is  derived  from 
the  data  of  Klstler^1^  Instead  of  by  a  mixing  length  concept.  The  normal 
pressure  gradient  in  the  boundary  layer  due  to  turbulence  Is  taken  into 
account  by  Including  the  y  momentum  equation  of  mean  motion.  The  fluctuation 
term  <(pv)‘v'>  In  the  y  momentum  equation  is  correlated  with  the  Reynolds  shear 
stress  to  predict  the  pressure  variation  within  the  boundary  layer.  The 
sensitivity  of  the  temperature  distribution  and  heat  transfer  due  to  the 
Influence  of  a  turbulent  Prandtl  number  profile  variation  Is  examined.  Pre¬ 
dictions  of  the  velocity,  static  temperature  profiles; as  well  as  the  wall  skin 
friction  and  heat  transfer  values,  are  compared  with  experimental  data. 
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SECTION  II 


GOVERNING  EQUATIONS 


The  basic  equations  are  written  In  essentially  the  form  as  given  by 
Van  Drlest^1^  except  for  the  additional  normal  momentum  equation 
Continuity 

^  +  ^  [pv  +  <(pv)’>]  «  0 
Streamwlse  momentum  equation 

pU  g  ♦  [PV  *  «<■»)•>]  |  ■  -  S  ♦  %  <*fc>  •  #  [<(PV)'U'>] 

Normal  momentum  equation 

=  0 


pu 


jg  [<(pv)V>  +p] 


Energy  equation 

&  ♦  [p, ♦  <■»•>)  §  *  #  [f  | 


3u 


(1) 


(2) 


(3) 


(4) 


To  perform  calculations  of  turbulent  boundary  layers.  Initial  conditions  are 

required  together  with  the  usual  nonsllp  condition  and  matching  temperature 

condition  at  the  wall.  The  boundary  conditions  at  the  outer  edge  of  the 

boundary  layer  have  u  and  h  approaching  ufi  and  hg,  respectively. 

(121 

The  fluctuation  terms  can  be  expanded  as  follows:1  7 

< (pv ) ' u ' >  =  p<u'v*>  ♦  v<p'u'>  +<j/u'v'>  (5) 

< ( pv ) ' v ' >  ~  p<v'v*>  +  v<p'v'>  +<p'v'v'>  (6) 

< ( pv ) *  h ' >  =  p<v'h'>  ♦  v<p'h’>  +<p'v'h'>  (7) 

The  above  turbulent  correlation  terms  are  not  known  a  priori.  If  the  above 
equations  are  normalized  with  respect  to  a  reference  condition,  an  estimation 
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of  the  aforementioned  terms  can  be  made.  For  a  boundary  layer  this  may  be 
done  with  a  standard  order  of  magnitude  analysis.  Herring  and  Mellor^) 
have  shown  the  order  of  magnitude  for  these  turbulent  correlations.  In  their 
analysis,  Ap  and  Ah  denote  the  variation  of  mean  density  and  mean  static 
enthalpy,  respectively,  across  the  boundary  layer.  Ap/pfi  is  restricted  to  the 
order  of  magnitude  of  unity.  The  second  term  contains  v  In  Eqs.  (5),  (6)  and 
(7),  and  v/ue  Is  proportional  to  6/L.  Therefore,  v<p'u'>/ue2Ap  has  the  order 
of  (6/L)2,  where  6  Is  the  boundary  layer  thickness  and  L  Is  the  length  of  the 
plate.  These  terms  are  neglected  by  Van  Driest^2)  for  a  thin  boundary 
layer.  The  order  of  magnitude  of  the  triple  correlation  terms  in  Eqs.  (5), 

(6)  and  (7)  Is  difficult  to  assess.  For  the  present  purpose,  an  estimate 

by  their  root-mean- square  values  seems  to  be  adequate.  We  have 
<p Vv'>/pue2-.0(/“pr?6/pL).  A  similar  argument  leads  to  the  Identical 
results  for  Eqs.  (6)  and  (7).  The  upper  limit  of  /p'"7p  can  be  obtained  from 
the  data  of  Kistler  (p'/p  *  -  T'/T).  The  relative  density  fluctuation  Is 
around  one  tenth  of  2(TQ-  TJ/(T0  +  Tw)  for  Mach  numbers  up  to  around  five. 

In  this  form  the  relative  density  fluctuation  is  nearly  independent  of  Mach 
number.  In  view  of  the  order  of  magnitude  analysis,  the  only  triple  correla¬ 
tions  retained  In  the  present  analysis  are  <p*u'v'>,  <p'vV>  and  <pVh'>  . 
Viscosity  fluctuation  Is  neglected  following  the  work  of  Bushnell  et  al.^ 

The  Reynolds  shear  stress  becomes 

<(pv)'u‘>  *  p<u'v’>  ( 1  +apuv  )  (8) 

Similar  results  can  h  expressed  for  Eqs.  (6)  and  (7) 


<(pv)  V>  =  p<v'v‘>  (l  +apvv 

(9) 

<(pv) ' h * >  =  p<v'h'>  (l  +apyh 

(10) 
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where  the  a's  are  the  correlation  functions.  The  subscripts  denote  the 
respective  variables  to  be  correlated.  Since  the  density  fluctuation  is 
introduced  by  phenomenlogical  means  with  uncertainty  in  the  RMS  value,  the 
correlation  coefficients  are  assumed  to  be  unity.  This  assumption  has  been 
applied  implicitly  in  the  work  of  Bushnell  et  al.^  for  hypersonic  turbulent 
boundary  layers. 

In  the  framework  of  the  eddy  viscosity  concept,  the  apparent  viscosity 
coefficient  becomes 


-  p<u'v‘>  -  <p'uV>  =  pe(l  + 


(ID 


The  turbulent  Prandtl  number  can  be  defined  as 


Pr, 


<u 

<v 


v>('  *  %,  )  (%) 
'T,>(' ♦  Ch  p  )  (§ ) 


(12) 


In  the  present  analysis,  the  turbulent  Prandtl  number  reduces  to  the 


conventional  form 


(4,13) 


by  assigning  opuy  »c.pyh  *  1. 


The  apparent  mass  flux  <p'v*>  appears  consistently  with  pv  in  the  same 
form  throughout  the  governing  equations.  This  term  is  easily  taken  into 
account  b«  defining 


v  * 


V  + 


<p'v'> 

P 


(13) 


The  correlation  term  p<v'v'>  in  the  normal  momentum  equation  is  closely  related 
to  the  turbulent  shear  stress.  Bradshaw  has  shown  that  the  turbulent 
shear  stress  can  be  approximated  as  follows: 

<uV>  =  0.15  (u12  +  v'2  +  w'2)  (14) 

(151 

The  correlation  can  be  observed  from  the  data  of  Klebanoff.  '  A 
similar  correlation  can  be  observed  between  <v,2>  and  the  turbulent  kinetic 
energy,  (u'2  +  v'2  +  w'2).  According  to  the  data  of  Ref.  15,  the  correlation 


5 


constant  appears  to  be  about  0.2  over  the  major  portion  of  the  boundary  layer. 
The  profiles  of  <u'v’>  /q  and  <v'2>/q  across  a  turbulent  boundary  layer  can 
be  deduced  from  the  data  of  Klebanoff^)  and  is  presented  in  Figure  1. 

<V2>=  0.2  (u'2  +  v12  +  w'2)  (15) 

Therefore,  the  correlation  term  of  p<v'v'>  may  be  approximated  as 

p<vV>  =1.4  pe  (16) 

Finally,  closure  of  the  governing  equations  was  obtained  by  introducing  the 
perfect  gas  law  and  Sutherland's  viscosity  equation.  The  Levy-Lees 
transformation  converts  the  governing  equations  into 


vn  ♦  2CF?  ♦  F  .  0 

(17) 

(CcF  )  -  vF  +  B(e-F’)  . 
n  n  n 

2CFF^ 

(18) 

r  „  cu  (e-  d  ■ 

L  *  Pe  (2CP 

=  0 
n 

(19) 

(pFce0n)  ■*,'  c?F*n 
'n 

-  2^F? 

(20) 

where  v  is  defined  as 


V 


IL 


p  U  u 

Ke  e  e 


The  set  of  equations  is  identical  to  that  of  Ref.  5,  except  for  the  additional 
normal  momentum  equation,  which,  however,  indicates  that  the  variation  of  the 
static  pressure  across  the  turbulent  boundary  layer  is  small.  The  pressure 
variation  is  such  that  p  +  p<v'2>(l  +  /p‘ J/p)  remains  constant  within  the 
turbulent  boundary  layer.  At  the  wall  <v'2>  vanishes  and  there  the  surface 
pressure  then  reaches  the  maximum  value. 

The  numerical  scheme  used  in  the  present  analysis  was  identical  to  that  of 
Ref.  5,  except  the  density  variation  was  deduced  from  the  equation  of  state  by 
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the  computed  temperature  and  pressure.  The  continuity  equation  was  solved  by 
a  straightforward  differencing  scheme.  The  normal  momentum  equtlon  reduced 
to  a  single  algebraic  equation  related  to  the  local  Reynolds  shear  stress. 

The  streamwi^e  momentum  equation  and  energy  equation  were  solved  simultaneously 
by  an  Implicit  scheme  together  with  an  Iteration  procedure  of  the  nondlmenslonal 
velocity  and  temperature  profiles.  The  detailed  description  of  the  numerical 
method  Is  Included  In  Refs.  5  and  16. 


SECTION  III 
EDDY  VISCOSITY  MODEL 

Several  multilayer  viscosity  models  have  been  used  successfully  to 
predict  the  flow  properties  of  compressible  turbulent  boundary  layers. 

The  differences  among  numerical  solutions  by  various  eddy  viscosity  models 
are  small. ^  The  crucial  po  nt  In  selecting  the  eddy  viscosity  model  Is  the 
damping  factor  In  the  viscous  sublayer,  which  determines  the  length  scale,  skin 
friction  and  heat  transfer.  The  most  versatile  damping  factor  in  which 
the  effect  of  mass  transfer  can  be  easily  included  by  adjusting  a  constant, 
is  due  to  Van  Driest.  In  fact,  Bushnell  et  al.^  have  established  the 
correlation  between  the  damping  constant  and  the  dimensionless  blowing  rate. 
The  two-layer  eddy  viscosity  model  for  the  hypersonic  turbulent  boundary  layer 
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can  be  given  as 


Inner  region 


(21) 


where  kj  is;  the  well-known  von  Karman  constant.  In  the  present  analysis,  a 
value  of  0.4  was  used  for  kx.  The  inner  layer  model  differs  from  that  of 
Cebeci's  work'  '  only  in  the  evaluation  of  the  shear  stress  in  the  damping  factor. 
In  the  present  analysis  the  shear  stress  assumes  the  local  value  instead  of  the 
wall  value  used  in  Cebeci's  work.  The  last  term  in  Eq.  (21) is  the  modification 
for  density  fluctuation  at  hypersonic  Mach  numbers.  The  density  fluctuation 
term  also  appears  in  the  outer  region  of  the  eddy  viscosity  model.  We  have 

ft)  ■  (m,«V(J)/»)  («  ♦f1)  <22> 

o 

The  constant  k2  is  given  the  value  of  0.0168.  The  quantity  y(y/6)  is  the  so 

f6 

called  intermittency  factor  in  the  law  of  thewake  region.  The  6?  (l-u/ug)dy 
is  the  conmonly  accepted  length  scale  from  the  velocity  defect  law. 

The  transition  model  of  Dhawan  and  Harasimha^^  was  adopted  in  the 
present  analysis 


T  ‘P(v  ♦  re)  !£ 

where 

r(x)  *  1  -  exp  (-.412  x2) 

The  normalized  streamwise  coordinate  in  the  transition  region  is  defined  as 


The  parameter  A  is  a  measure  of  the  extent  of  the  transition  region  defined  by 

A  =  xr«3A  "  xr=l/„ 
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The  transition  model  successfully  described  the  transition  behavior  for  the 

1 5) 

supersonic  compressible  turbulent  boundary  layer.'  '  The  transition  model 
requires  knowledge  of  the  beginning  of  the  transition  and  the  extent  of  the 
transition  region.  According  to  a  substantial  amount  of  experimental  data, 
the  Reynolds  number  at  the  end  of  the  transition  region  is  t, ice  its  value  at 
the  beginning  of  transition.  This  experimental  observation  was  Incorporated 
into  the  transition  model.  Therefore  the  only  information  required  was  the 
onset  of  the  location  of  transition.  For  all  the  calculated  cases,  this 

•j  ■> 

location  is  contained  by  1.75  x  10  <  ReQ  <  2.2  x  10  . 

—  U  — 

The  present  modification  of  the  eddy  viscosity  model  for  extension  to 
hypersonic  flows  is  accomplished  by  including  the  density  fluctuation  term. 
Bushnell  et  al.  introduced  the  density  fluctuation  term  by  a  mixing  length 
type  model.  Essential  agreement  between  their  numerical  solution  and  hyper¬ 
sonic  turbulent  boundary  layer  data  was  obtained.  Their  achievement  confirms 
the  basic  contention  that  the  density  fluctuation  term  at  hypersonic  speeds  is 
significant.  The  density  fluctuation  term  of  the  present  analysis  was 
obtained  from  the  data  of  Kistler  at  Mach  numbers  up  to  4.7  and  for  near 
adiabatic  wall  conditions.  This  modification  was  then  used  to  perform  cal¬ 
culations  at  higher  Mach  numbers  and  with  heat  tranrfer.  The  justification 
is  based  on  the  fact  that  the  relative  temperature  fluctuation  scaled  by 
2(Tq-  TJ/(T  +  T^)  is  nearly  independent  of  Mach  number.  Kovansnay 
Kistler^ ^  and  Laderman  et  al.^  have  found  a  strong  negative  correlation 
between  the  temperature  and  the  velocity  fluctuation  <u'T’>.  The  static 
temperature  fluctuation  levels  are  proportional  to  the  temperature  difference 
across  the  boundary  layer. Further,  the  temperature  fluctuation  level 
obtained  by  Laderman  and  Demetriades  for  a  highly  cooled  wall  at  a  Mach  number 
of  9.37  was  similar  to  the  dat^  of  Kistler,  although  at  a  lower  value. 
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From  ihe  equation  of  state,  the  fluctuating  thermodynamic  properties  :an 
be  expressed  as 


The  temperature  fluctuations  obtained  by  Klstler^1^  were  based  on  the 
assumption  that  the  pressure  fluctuation  Is  negligible.  Thus  the  temperature 
and  density  fluctuations  are  identical  in  magnitude  but  opposite  In  sign.  This 
equality,  which  also  has  been  pointed  out  by  Laufer.09)  jS  not  a  direct  result 
of  measurement. 

In  the  present  analysis,  a  pressure  fluctuation  term  has  been  obtained 
from  the  normal  momentum  Eq.  (3).  The  origin  of  this  pressure  variation  is 
attributed  by  the  normal  component  of  the  Reynolds  shear  stress.  The  corre¬ 
lated  pressure  variation  is  proportional  to  yMe2(2x/pu2).  It  is  interesting 
to  note  that  the  prerent  result  is  almost  identical  to  the  experimental  result 

of  Kistler  et  al.'a^  They  estimated  the  relative  level  of  the  pressure 
fluctuation  to  be  <p>/p  =2.5  Mg2Cf  for  Mg  >  2.  Further,  Laderman  and 
Demetr1ade^22^  have  shown  that  the  relative  magnitude  of  the  density  fluctua¬ 
tion  is  greater  than  that  of  the  pressure  and  temperature  fluctuation.  In  order 
to  be  consistent  with  the  present  formulation,  the  RMS  value  of  the  density 
fluctuation  is  tentatively  given  as 

@L  ,  +  .<(P.vr..y.,>  (23) 

P  T  p 

The  RMS  value  of  the  temperature  fluctuation  is  introduced  as 

*T 1  *  2  (fc£)  f  (f)  <24> 

where  f(y/5)  gives  a  fitting  of  Kistler's  data,  having  a  maximum  value  of 
0.1.  The  pressure  fluctuation  term  was  calculated  from  the  normal  momentum 
equation. 
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SECTION  IV 


TURBULENT  PRANDTL  NUMBER 

Recently  Simpson  et  al.^®^  have  established  the  uncertainty  envelope 

of  the  turbulent  Prandtl  number  for  Incompressible  zero  pressure  gradient 

turbulent  boundaries.  They  conclude  that  in  the  wall  region  the  molecular 

viscosity,  Prandtl  number  and  the  small  scale  turbulence  govern  the  momentum 

and  heat  transport.  Furthermore,  the  turbulent  Prandtl  number  Is  greater  than 

unity  at  the  wall.  However,  In  the  outer  region  Prt  is  less  than  unity. 

Several  turbulent  Prandtl  number  models  have  been  developed. The  most 

common  characteristic  is  that  the  predicted  Pr^  value  exceeds  unity  near  the 

wall.  The  maximum  value  of  Pr^  seems  to  be  1.35  at  the  wall.  Experimental 

data  by  Meier  and  Rotta  at  Mach  numbers  up  to  4.5  and  near  adiabatic  wall 

conditions  yield  a  Pr^  distribution  with  a  value  greater  than  unity  near  the 

wall.  Data  by  Horstman  and  Owen^20^  at  Me  *  7.2  and  cooled  wall  conditions 

also  indicate  this  behavior.  Nevertheless,  all  the  predictions  fall  into  the 

uncertainty  envelope  of  Simpson  et  al.  In  view  of  the  unresolved  question  on 

(21 1 

the  Prt,  an  empirical  Prt  distribution  has  been  suggested  by  Rotta'  '  as 
follows: 

Prt  =  .95  -  0.45  (y/$)2 

This  appears  to  be  a  very  good  approximation  of  the  Prt  in  the 
outer  region.  For  the  purpose  of  studying  the  sensitivity  of  Pr^  on  the 
present  numerical  solutions,  a  similar  empirical  formulation  is  suggested  as 
follows: 

Prt  =  Pr,  e"10(y/6)  +  Pr2  [l.  -  0.2  (£)]  (25) 

where  0.8  <  Pr2  <  1.0,  .2  <  Pr!  <  .4. 
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The  above  approximation  closely  follows  the  outer  bounds  of  the  data. 
If  Pr2  Is  assigned  to  be  0.90,  Prx  s  0.30,  the  Pr^  distribution  by  Eq.  (25) 
describes  approximately  the  mean  value  of  the  uncertainty  envelope.  A 
graphical  presentation  of  the  uncertainty  envelope  of  the  turbulent  Prandtl 
number  and  the  present  approximations  of  its  limits  are  given  in  Figure  2. 


SECTION  V 

DISCUSSION  OF  RESULTS 

The  computed  results  are  presented  In  three  sections.  The  first  portion 
of  the  results  is  concerned  with  the  pressure  variation  due  to  turbulence 
across  the  hypersonic  turbulent  boundary  layer  and  Its  influence  on  the  char¬ 
acteristics  of  the  shear  layer.  The  second  section  of  the  presentation  is  an 
evaluation  of  the  turbulent  heat  flux  In  terms  of  the  Prandtl  number  variation. 
This  parameter  Is  a  measure  of  the  relative  magnitude  between  the  eddy 
viscosity  and  the  apparent  heat  conductivity  Pr^  =  Cppe/At.  The  remaining 
section  Is  devoted  to  the  solution  of  hypersonic  turbulent  boundary  layers 
with  the  inclusion  of  the  density  fluctuation  terms.  The  validity  of  this 
eddy  viscosity  model  In  the  hypersonic  flow  regime  with  heat  transfer  Is 
determined  by  means  of  comparison  with  experimental  data. 
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The  third  figure  presents  the  static  pressure  variation  across  the 
turbulent  boundary  layer  at  Mach  number  9.37.  The  experimental^  data  were 
obtained  by  Laderman  and  Oemetriades. ^  The  normalizing  boundaryVlayer 
thickness  in  the  presentation  of  their  data  Is  considered  to  be  4.4  inches 
according  to  their  revised  value. ^22^  The  calculated  pressure  distribution 
within  the  boundary  layer  agrees  very  well  with  the  measurement.  The  pressure 
decreases  very  rapidly  from  the  wall,  reaches  the  minimum  value  near  the  outer 
edge  of  the  sublayer,  and  then  approaches  the  free  stream  value  asymptotically. 
The  pressure  variation  across  the  turbulent  boundary  is  directly  proportional 
to  the  turbulent  shear  stress  distribution.  The  relative  magnitude  of  the 
pressure  variation  is  proportionally  to  the  square  of  the  free  stream  Mach 
number  and  the  correlation  function  between  <v'2>  and  the  turbulent  kinetic 
energy.  The  present  modeling  of  the  pressure  variation  assumed  no  free 
stream  turbulence.  Therefore  the  surface  pressure  equals  the  free  stream 
value.  If  the  free  stream  turbulence  value  can  be  assessed,  the  maximum 
pressure  will  be  attained  at  the  surface  as  that  reported  by  Fischer  et  al.v  ' 
In  order  to  evaluate  the  significance  of  the  small  pressure  variation 
on  the  mean  motion,  calculations  were  performed  with  and  without  the  inclusion 
of  the  pressure  fluctuation  term.  The  difference  between  the  two  calculations 
is  negligible  for  the  cases  investigated.  In  addition,  a  reduction  in 
streamwlse  step  size  for  the  numerical  scheme  must  be  Implemented  to  Include 
the  normal  momentum  equation.  In  spite  of  the  small  transverse  pressure 
variation,  the  numerical  procedure  required  a  significantly  large  number  of 
iterations  to  meet  the  established  convergent  criterion.  The  solution  even 
failed  to  converge  if  too  coarse  a  streamwise  mesh  size  were  assigned.  Based 
upon  the  aforementioned  observations,  the  small  pressure  variation  across  the 
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turbulent  boundary  layer  was  neglected  In  the  mean  equation  of  motion  In  the 
following  analysis. 

To  examine  the  Influence  of  turbulent  Prandtl  number  on  numerical 
solutions,  an  empirical  formulation  of  the  turbulent  Prandtl  number  was  int  o • 
duced.  The  empirical  expression  provided  a  close  description  of  the  upper 
and  lower  limits  of  the  uncertainty  envelope  of  the  turbulent  Prandtl  number. 
The  molecular  Prandtl  number  was  assigned  the  value  of  0.73  throughout  the 
test.  Holden's  data^)  of  the  skin  friction  coefficient  and  Stanton  number 
were  used  for  the  evaluation  process.  His  turbulent  flow  data  were  taken 
under  a  natural  laminar- turbulent  transition  condition. 

Figure  4  presents  the  skin  friction  coefficient  and  Stanton  number 
distributions  at  *  7.97  with  a  wall  to  free  stream  stagnation  temperature 
ratio  of  T^/T  *  0.308.  Two  numerical  solutions  are  presented  In  the  graph. 
One  of  the  solutions  Is  computed  by  assuming  the  upper  11m<t  of  the  Prandtl 
number  variation,  and  the  other  assumes  the  lower  limit.  In  the  present 
context,  the  Pr^  distributions  across  the  boundary  layer  were  generated 
from  Eq.  (25).  The  upper  limit  of  the  Prt  variation  Is  obtained  by  assigning 
Pr2  *  0.4  and  Pr2  *  1.0;  the  lower  limit  of  the  Prt  variation  Is  calculated 
with  Prj  *  0.2  and  Pr2  *  0.8.  The  difference  In  the  boundary  layer  calcula¬ 
tions  was  found  to  be  small.  The  computed  Stanton  number  and  skin  friction 
coefficient  distributions  indicate  excellent  agreement  with  the  experimental 
data.  The  difference  between  the  two  solutions  is  noticeable  downstream  of 
the  transition  region.  The  solution  associated  with  the  upper  limit  Pr^ 
variation  yields  a  lower  value  of  than  the  solution  by  the  lower  limit  Prt> 
However,  the  trend  Is  reversed  In  predicting  the  Stanton  number.  As  a 
consequence,  the  Reynolds  analogy  factor  2St/Cf  Is  decreased  with  decreasing 
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turbulent  Prandtl  number.  The  Identical  observation  can  be  made  In  the  the 
work  of  Bushnell  and  Beckwith. 

The  Influence  of  turbulent  Prandtl  number  on  the  detailed  flow  field  was 
also  Investigated  for  the  two  limits  of  Prt  distributions.  The  computed 
velocity  profiles  are  nearly  Identical,  thus  only  static  temperature  profiles 
are  presented  In  Figure  5.  The  comparison  between  the  data  of  Laderman 
et  al.^  and  the  computed  results  Is  fair.  One  observes  that  the  present 
calculations  predict  correctly  the  location  of  the  maximum  temperature  In  the 
turbulent  boundary  layer.  The  p<*ak  temperature  Is  located  near  the  outer  edge 
of  the  viscous  sublayer.  The  difference  between  the  two  solutions  of  the 
limiting  Prandtl  number  distributions  Is  small.  Significant  differences 
between  the  two  solutions  appear  only  In  the  magnitude  of  the  maximum 
temperature  and  In  the  law  of  the  wake  region. 

Calculations  with  a  constant  Pr^  value  of  0.9  and  the  Prt  distribution 
suggested  by  Rotta^21^  have  also  been  performed.  The  calculations  are  con¬ 
tained  within  the  limiting  solutions  In  either  C^,  St  or  velocity  and  static 
temperature  profiles.  In  particular,  the  small  difference  (less  than  3%) 
between  solutions  of  Prt  ■  0.9  and  Rotta's  Pr^  distribution  Indicates  an 
Insensitivity  of  the  numerical  solutions  to  the  Pr^  variation  In  the  law  of  the 
wake  region. 

In  order  to  ascertain  the  predicted  trend  of  the  skin  friction  and  the 
Stanton  number  due  to  different  values  or  Pr^,  two  additional  calculations 
with  constant  turbulent  Prandtl  numbers  of  0.7  and  1.3  at  Me  *  10.57  have 
been  Included  In  the  present  study.  In  Fig.  6,  the  calculated  velocity  and 
temperature  profiles  at  different  but  constant  turbulent  Prandtl  numbers  are 
presented.  One  observes  that  the  velocity  profile  Is  hardly  affected  by  the 
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different  magnitude  of  the  constant  turbulent  Prandtl  number.  On  the  other 
hand,  significant  disparity  In  temperature  profiles  Is  obvious.  The  maximum 
deviation  between  two  solutions  appears  in  the  magnitude  of  the  predicted 
peak  temperature  and  the  temperature  profile  In  the  law  of  the  wake  region. 
This  behavior  Is  also  perceptible  In  Fig.  5  for  variable  Pr^  distributions  at 
different  values.  The  calculation  by  the  higher  Pr^  value  yields  a  higher 
static  temperature  than  the  calculation  with  a  lower  value  of  Prt  in  the 
laminar  sublayer  region.  As  a  consequence,  the  calculation  by  the  higher 
value  of  turbulent  Prandtl  number  generates  a  lower  eddy  viscosity  for  the 
entire  Inner  region  (Figure  7).  Therefore,  the  solution  with  Pr^  *  1.3  yields 
a  higher  heat  transfer  rate  but  lower  shear  stress  at  the  wall  than  the 
solution  with  Pr^  *  0.7.  In  all,  the  Stanton  number  and  the  skin  friction 
coefficient  exhibit  a  relative  Insensitivity  with  respect  to  the  turbulent 
Prandtl  number  variation.  The  difference  of  46X  In  Pr^  produces  a  difference 
of  only  7%  in  the  Stanton  number  and  4%  in  C^. 

The  conventional  comparison  of  (H-Hw)/(He~Hw)  vs  u/ug  is  given  in  Fig.  8 

with  the  data  of  Laderman  and  Demetriades  for  *  9.37.  In  these  coordinates 

e 

Crocco's  relationship  for  Prt  *  Pr  ■  1  is  a  straight  line  which  Is  verified 
by  the  present  calculation.  The  lower  value  of  Prt  In  the  present  calculation 
produces  a  profile  closer  to  the  quadratic  law  than  the  higher  value  of  Pr^.. 
The  discrepancy  betweei  the  data  and  the  present  calculation  may  be  attributed 
to  the  aspect  of  flow  history.  The  deviation  from  the  Crocco  relationship 
has  been  documented  for  data^’^  obtained  from  wind  tunnel  walls.  Unfor¬ 
tunately,  detailed  upstream  flow  conditions  are  not  available,  and,  hence, 
it  is  not  possible  to  verify  this  contention. 
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In  Fig.  9»  the  evaluation  of  the  present  eddy  viscosity  model  with  an 

explicit  density  fluctuation  dependence  was  presented  together  with  the  data 

by  Holden  and  the  data  by  Hopkins  et  al .  The  experimental  data  were 

obtained  at  a  Mach  number  of  about  7.4  and  T '  /T  from  0.172  to  0.418.  All 

w  o 

data  were  obtained  under  natural  laminar  to  turbulent  transition  conditions. 
Holden's  data  of  the  Stanton  number  provided  a  distinctive  onset  location  of 
the  transition  zone.  Two  numerical  solutions  consistent  with  the  transition 
condition  were  obtained.  The  third  calculation  was  obtained  by  assigning 
transition  at  the  origin.  The  present  eddy  viscosity  model  which  Includes  the 
density  fluctuation  terms  predicts  very  well  the  skin  friction  coefficient. 

The  conventional  eddy  viscosity  model  without  the  density  fluctuation  correc¬ 
tion  underpredicts  Cf  in  the  low  Re0  region  but  Improves  steadily  downstream 
of  the  transition  zone.  The  calculated  Stanton  number  uniformly  overpredicts 
the  data  of  Holden.  The  difference  between  solutions  of  the  eddy  viscosity 
model  without  the  density  fluctuation  correction  and  the  present  model  Is  7%. 
The  turbulent  Prandtl  number  used  In  these  calculations  was  assumed  to  have  a 
constant  value  of  0.9.  This  simplification  was  made  In  the  present  calcula¬ 
tions  so  that  a  consistent  comparison  could  be  achieved  with  several  other 
investigations/  * 

Similar  comparisons  at  Mach  number  around  eight  are  presented  In 
Figure  10.  Excellent  agreement  Is  Indicated  between  the  heat  transfer  data  of 
Holden  and  the  present  calculations  In  the  regions  of  laminar,  transition  and 
turbulent  flows.  The  eddy  viscosity  model  with  the  density  fluctuation  cor¬ 
rection  underpredicts  the  data  of  Hopkins  et  al.  In  the  low  Re.  region.  The 
eddy  viscosity  model  without  the  density  fluctuation  correction  would 
underpredict  the  data  further  by  8%.  Agreement  between  the  data  and  the 
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calculated  Cf  by  the  eddy  viscosity  model  with  the  density  fluctuation 
correction  Is  nearly  attained  In  the  higher  Re0  region.  The  solution  by  the 
eddy  viscosity  model  without  the  density  fluctuation  correction  slight 
underpredicts  the  Cf  and  St  measured  values  In  the  higher  ReQ  region. 

In  Fig.  11,  the  calculations  are  compared  with  the  data  of  Holden  at  a 
free  stream  Mach  number  around  ten  and  TW/TQ  about  0.2.  The  numerical  solution 
by  the  eddy  viscosity  model  without  the  density  fluctuation  correction  under¬ 
predicts  the  skin  friction  coefficient  over  the  complete  range  of  Re..  These 

y 

calculated  heat  transfer  rates  show  better  comparison  with  the  measured  results 
than  that  by  the  eddy  viscosity  model  with  the  density  fluctuation  correction. 
However,  the  prediction  of  the  heat  transfer  in  the  laminar  and  transitional 
regions  Is  excellent.  The  eddy  viscosity  model  with  the  density  fluctuation 
correction  yields  essential  agreement  with  the  Cf  measurements  over  the 
entire  range  of  Re0.  An  Interesting  observation  may  be  noted  In  Fig.  11; 
i.e.,  turbulent  toundary  layers  at  low  Reynolds  numbers,  exhibit  relatively 
high  values  of  the  shear  stress  and  heat  transfer  at  the  wall.  This  phenomena 
can  be  detected  for  all  the  cases  Investigated.  The  behavior  Is  particularly 
noticeable  by  comparing  the  data  and  the  solution  with  an  eddy  viscosity 
model  with  the  density  fluctuation  correction.  The  numerical  solution  always 
underpredicts  Cf  and  St  for  the  lower  ReQ  region  but  reaches  near  agreement 
at  the  higher  ReQ  region.  The  phenomenon  Is  progressively  more  pronounced  as 
the  free  stream  Mach  number  increases. 

In  Fig.  12,  the  calculations  by  an  eddy  viscosity  model  with  and  without 
the  density  fluctuation  correction  are  compared  with  the  data  by  Holden. 

His  data  were  obtained  at  =  12.04  and  TW/TQ  *  0.157.  The  maximum  value 

of  ReQ  for  the  experiment  Is  6.849  x  lO'*.  The  prediction  In  heat  transfer 

0 
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Is  again  excellent  In  the  laminar  and  transition  regions.  Near  agreement 
Is  also  attained  in  the  turbulent  region.  The  skin  friction  coefficient 
calculation  is  also  underpredicted  In  this  relatively  low  Re.  region.  The 

U 

calculation  with  the  eddy  viscosity  model  without  the  density  fluctuation 
correction  further  underpredicts  by  10%.  To  substantiate  the  present  results, 
those  of  Van  Driest's  theory  are  also  presented  in  Figure  12.  The  difference 

between  the  predictions  of  Van  Driest  theory  and  the  present  model  with  the 
density  fluctuation  correction  is  small.  The  measured  data  indicate  a 

higher  value  than  both  calculations  In  the  low  Re0  region. 

A  summary  of  the  eddy  viscosity  model  with  the  density  fluctuation 

correction  over  a  wide  range  of  the  free  stream  Mach  number  Is  presented  in 

(21 

Figure  13.  Van  Driest's  theory  suggested  by  Hopkins  et  al.  #  is  adopted  to 
be  used  as  the  criterion.  The  comparison  is  divided  into  two  groups.  The 

4 

first  group  is  restricted  to  Re0  less  than  10  ,  while  the  second  group  of 
comparisons  is  presented  for  Re0  greater  than  10*.  The  division  of  Re0  in  the 
present  comparison  is  rather  arbitrary.  The  noticeable  underprediction  of 
Cf  in  the  low  Re0  region  is  obvious.  The  difference  between  solutions  by  the 
eddy  viscosity  model  without  the  density  f  ijctuation  correction  and  the  pre¬ 
sent  model  with  density  fluctuation  correction  diminishes  as  the  free  stream 
Mach  number  decreases.  The  difference  becomes  negligible  as  the  Mach  number 
reaches  a  value  less  than  two.  This  behavior  is  expected  because  the  density 
fluctuation  correction  Is  proportional  to  TQ-Te.  The  difference  between  the 
Van  Driest  theory  and  the  present  calculation  is  small,  (In  general,  both 
underpredict  in  the  lower  Re0  region)  but  is  within  the  scattering  of  the 
measurements  in  the  higher  Re0  region. 
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The  particular  behavior  of  the  turbulent  boundary  layer  at  low  Reynolds 
numbers  has  been  Investigated  by  Huffman  and  Bradshaw.  They  suggested 
that  a  turbulent-irrotatlonal  interface  exerts  a  significant  effect  on  the 
outer  region  of  a  boundary  layer  at  low  Reynolds  numbers.  They  also  indicated 
that  one  of  the  primary  effects  of  external  Influences  on  the  inner  layer  is  a 
change  in  the  damping  constant  A.  In  Fig.  14,  the  dependence  of  the  skin 
friction  coefficient  on  the  damping  constant  A  is  presented.  The  original 
value  of  A  in  the  Van  Driest  damping  factor  is  26.  A  value  range  of  A  from 
0.01  to  99  is  computed.  In  the  present  numerical  procedure,  the  smaller  value 
of  A  is  equivalent  to  letting  the  damping  factor  D  approach  unity.  Thus,  the 
law  of  the  wall  is  extended  closer  to  the  surface.  On  the  other  hand,  the 
greater  value  of  A  yields  a  value  of  which  approaches  that  of  the  laminar 
flow.  In  this  given  value  range  of  A,  the  calculated  varies  by  a  factor  of 
three.  The  deficiency  in  predicting  C-  at  the  lower  Re^  region  can  be  easily 
corrected  by  either  adjusting  the  sublayer  thickness  through  the  damping 
constant  or  by  removing  the  i ntermi ttency  correction  in  the  law  of  the  wake 
region.  However,  the  understanding  of  the  phenomenon  can  be  obtained  only 
by  extensive  experimental  investigations. 
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SECTION  VI 


CONCLUDING  REMARKS 

Numerical  solutions  to  the  mean  equations  of  turbulent  motion  were 
obtained  by  a  finite-difference  scheme.  Models  of  the  terms  for  the  Reynolds 
stress  and  the  turbulent  enemy  flux  were  modified  to  Include  the  effect  of 
fluctuating  density  which  becomes  significant  at  hypersonic  speeds.  The 
Influence  of  the  turbulent  Prandtl  number  on  a  numerical  solution  was  also 
Investigated. 

The  present  eddy  viscosity  model  with  an  explicit  dependence  on  the 
density  fluctuation  term  predicts  the  skin  friction  coefficient  and  the  heat 
transfer  at  hypersonic  Mach  numbers.  The  agreement  between  data  and  the 
calculations  Is  excellent  in  the  laminar  and  transition  regions.  In  the 
turbulent  flow  regime,  the  present  solutions  as  well  as  all  current  theories 
underpredict  the  skin  friction  data  in  the  lower  Re0  region.  However,  the 
accuracy  Improves  steadily  as  the  flow  proceeds  farther  downstream  of  the 
transition  region.  In  the  higher  Re.  region,  the  agreement  between  the  data 
and  the  present  solutions  is  excellent.  On  the  other  hand,  the  eddy  viscosity 
model  without  the  density  fluctuation  correction  significantly  underpredicts 
the  skin  friction  coefficient  data.  The  discrepancy  between  data  and  calcula¬ 
tions  without  the  fluctuating  density  correction  becomes  progressively  more 
pronounced  as  the  free  stream  Mach  number  increases.  The  density  fluctuation 
correction  term  seems  to  be  nece*:st-ry  to  Improve  the  accuracy  of  solutions  In 
the  investigated  hypersonic  flow  regime. 

The  present  analysis  Indicates  that  in  a  turbulent  boundary  layer  the 
sum  p  +  <(pv)'v'>  remains  constant  across  the  shear  layer.  Therefore  the 
static  pressure  across  the  boundary  layer  must  adjust  Itself  to  accommodate 
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the  fluctuating  transverse  velocity  term.  According  to  the  present  model,  the 
maximum  magnitude  of  the  velocity  fluctuation  term  Is  about  five  percent  of 
the  free  stream  static  pressure  value  at  Mfi  =  9.37.  The  corresponding  static 
pressure  reaches  the  minimum  value  near  the  outer  region  of  the  laminar 
sublayer.  The  effect  of  this  pressure  variation  Is  negligible  for  the 
Investigated  cases,  however. 

The  numerical  solutions  exhibit  a  weak  dependence  on  the  turbulent 
Prandtl  number  variation.  For  turbulent  Prandtl  numbers  differing  by  more 
than  40%,  calculations  show  merely  a  6%  difference  In  and  Stanton  number. 
The  Reynolds  analogy  factor,  2St/C^  was  between  0.9  and  1.0  for  these  cases. 
These  encouraging  numerical  computations  show  that  the  density  fluctuation 
correction  Improves  the  overall  accuracy  of  the  solutions  In  the  hypersonic 
flow  region. 
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UNCERTAINTY  ENVELOPE 
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FIGURE  2  The  Uncertainty  Envelope  of  the  Turbulent  Prandtl  Number  and  the  Approximated  Prt  Distributions 
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FIGURE  3  Static  Pressure  Distribution  Across  Turbulent  Boundary  Layer 
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FIGURE  4  Effect  of  Variable  Turbulent  Prandtl  Number  on  C-  and  St 
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FIGURE  5  Effect  of  Variable  Turbulent  Prandtl  Number  on  Temperature  Profile 
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FIGURE  6  Effect  of  Constant  Turbulent  Prandtl  Number  on  Temperature  and  Velocity  Profil 
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FIGURE  7  Effect  of  Constant  Turbulent  Prandtl  Number  on  the  Eddy  Viscosity  Coefficient  Distribution 
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FIGURE  9  Effect  of  Density  Fluctuation  on  C*  and  St  at  M 
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FIGURE  10  Effect  of  Density  Fluctuation  on  Cf  and  St  a*  M  •*  8.0 
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FIGURE  11  Effect  of  Density  Fluctuation  on  C,  and  St  at  M  =  10.5 
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FIGURE  12  Effect  of  Density  Fluctuation  on  C-  and  St  at  M 


FIGURE  13  Skin  Friction  Coefficient  Variation  with  Mach  Number 
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FIGURE  14  Effect  of  Damping  Constant  on 
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